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We investigate the effects of geometrical micro-irregularities on the conversion efficiency of reactive 
flows in narrow channels of millimetric size. Three-dimensional simulations, based upon a Lattice- 
Boltzmann-Lajc-Wendroff code, indicate that periodic micro-barriers may have an appreciable effect 
on the effective reaction efficiency of the device. Once extrapolated to macroscopic scales, these 
effects can result in a sizeable increase of the overall reaction efficiency. 
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I. INTRODUCTION 

The formulation of mathematical models and atten- 
dant simulational tools for the description of complex 
phenomena involving multiple scales in space and time 
represents one of the outstanding frontiers of modern 
applied physics/mathematics 0]. One such example of 
complex multiscale phenomena is the dynamics of reac- 
tive flows, a subject of wide interdisciplinary concern in 
theoretical and applied science, with several applications 
in molecular engineering, material science, environmen- 
tal and life sciences alike. The complexity of reactive flow 
dynamics is parametrized by three dimensionless quanti- 
ties: the Reynolds number Re = UL/v, the Damkohler 
number Da — Th/Tc, and the Peclet number Pe — 
Here U , L and H denote the macroscopic flow speed and 
longitudinal/ transversal lengths of the flow, respectively, 
V the fluid kinematic viscosity and D the pollutant molec- 
ular diffusivity. The quantities Tc and th represent typical 
timescales of chemical and hydrodynamic phenomena. 

High Reynolds numbers are associated with turbu- 
lence, namely loss of coherence of the flow field in both 
space and time. High Damkohler numbers imply that 
chemistry is much faster than hydrodynamics, so that re- 
actions are always in chemical equilibrium and take place 
in tiny regions (thin flames, reaction pockets) of evolving 
flow configurations. The opposite regime ("well-stirred" 
reactor) characterizes situations where the chemistry is 
slow and always takes place at local mechanical equilib- 
rium. Finally, high Peclet numbers imply that the trans- 
ported species stick tightly to the fluid carrier (in the 
limit Pe — *■ oo the tracer field is "frozen-in" within flow 
streamlines). Navigation across the three dimensional 
Re — Da — Pe parameter space meets with an enormous 
variety of chemico-physical behaviours, ranging from tur- 
bulent combustion to hydrodynamic dispersion and oth- 
ers 1^. The picture gets further complicated when ge- 
ometry is taken into account, since boundary conditions 
select the spatio-temporal structures sustaining the non- 



linear interaction between the various fields. In this work 
we shall deal with low- Reynolds, fast-reacting flows with 
heterogeneus catalysis. In particular we wish to gain in- 
sights into the role of geometric micro-irregularities on 
the effective rate of absorption of tracer species (pollu- 
tant hereafter) at catalytic boundaries. This is a theme 
of broad interest, with applications in biology, physics, 
chemistry, environmental sciences and more. It is there- 
fore hoped that such kind of theoretical-computational 
studies may promote a better understanding of the com- 
plex phenomena behind these important applications [|j . 



II. MATHEMATICAL MODEL OF REACTIVE 
MICROFLOW DYNAMICS 

We shall deal with an incompressible, isothermal flow 
with soluted species which are transported (advect and 
diffuse) by the flow and, upon reaching solid walls, they 
undergo catalytic chemical reactions. The basic equa- 
tions of fluid motion are: 



dt{pu) 



dtp + div{pu) = (1) 
div{puu) = -VP + div{p\7u) (2) 



where p is the flow density, u the flow speed, P — pT 
the fluid pressure, T the temperature and p = pv the 
dynamic viscosity and uu denotes the dyadic tensor 

UaUh, a,b^ X, y, z. 

Multispecies transport with chemical reactions is de- 
scribed by a set of generalized continuity-diffusion equa- 
tions: 



dtCs + div(CsUs) = div{Ds^{Cslp)\ + tls 



(3) 



where Cg denotes the mass density of the generic s- 
th species, D^ its mass diffusivity and f^s is a surface- 
chemical reaction term to be detailed shortly. In the fol- 
lowing we indicate with the subscripts w and g the "wall" 
(solid) and "gas" in contact with the wall respectively. 
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According to Pick's law, the outgoing (bulk-to-wall) dif- 
fusive mass flux is given by: 



(4) 



Upon contact with solid walls, the transported species 
react according to the following empirical rate equation 
(the species index being removed for simplicity): 



n = 



dt 



(5) 



where the wall-flux is taken in the simple linear form: 

= if»(c»-c'") (6) 

where is the wall to/from fluid mass transfer rate 
and Kc is the chemical reaction rate dictating species 
consumption once a molecule is absorbed by the wall. 
The subscripts w and g mean "wall" (solid) and "gas" 
in a contact with the wall respectively. The above rate 
equation serves as a dynamic boundary condition for the 
species transport equations, so that each boundary cell 
can be regarded as a microscopic chemical reactor sus- 
tained by the mass inflow from the fluid. In the absence 
of surface chemical reactions the species concentration 
in the solid wall would pile up in time, up to the point 
where no outflow would occurr, a condition met when 
(jg _ (jw ^ Chemistry sets a time scale for this pile- up 
and fixes the steady-state mass exchange rate. At steady 
state we obtain: 



-CO 



hence 



C9 



(7) 



(8) 



where = l/K^ and Tc = l/Kc- These expres- 
sions show that finite-rate chemistry {Kc > 0) ensures 
a non-zero steady wall outflux of pollutant. At steady 
state, this mass flow to the catalytic wall comes into bal- 
ance with chemical reactions, thus fixing a relation be- 
tween the value of the wall-gradient concentration and 
its normal-to-wall gradient: 

where d± means the normal to the perimeter compo- 
nent of the gradient and p is the perimeter (volume / area) 
of the reactive cell. This is a mixed Neumann-Dirichlct 
boundary condition and identifies the free-slip length of 
the tracer as Ig = D{tw + Tc)Iv- 



III. THE COMPUTATIONAL METHOD 

The flow field is solved by a lattice Boltzmann method 
l^-Q while the multispecies transport and chemical re- 
actions are handled with a variant of the Lax-Wendroff 
method H . A few details are given in the following. 



A. Lattice Boltzmann equation 

The simplest, and most popular form of lattice Boltz- 
mann equation (Lattice BGK, for Bahtnagar, Gross, 
Krook) 0, reads as follows: 

/,(f -I- i + 1) - /,(x , t) = ~u{f, - ft\{x, t) (9) 

where fi{x,t) = f{x,v = Ci,t) is a discrete population 
moving along the discrete speed q. The set of discrete 
speeds must be chosen in such a way as to guarantee 
mass, momentum and energy conservation, as well as ro- 
tational invariance. Only a limited subclass of lattices 
qualifies. In the sequel, we shall refer to the nineteen- 
speed lattice consisting of zero-speed, speed one c — 1 
(nearest neighbor connection) , and speed c = ^/2, (next- 
nearest- neighbor connection). This makes a total of 19 
discrete speeds, 6 neighbors, 12 nearest-neighbors and 1 
rest particle (c — O). The right hand side of (^) represents 
the relaxation to a local equilibrium /f in a time lapse of 
the order of lu^^. This local equilibrium is usually taken 
in the form of a quadratic expansion of a Maxwellian: 



f! = P 
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II) 



2ci 



(10) 



where Cs is the sound speed and / denotes the identity. 

Once the discrete populations are known, fluid density 
and speed are obtained by (weighted) sums over the set 
of discrete speeds: 



pu 



(11) 



LBE was historically derived as the one-body kinetic 
equation resulting from many-body Lattice Gas Au- 
tomata, but it can mathematically obtained by standard 
projection upon Hermite polynomials of the continuum 
BGK equation and subsequent evaluation of the kinetic 
moment by Gaussian quadrature ||] . It so happens that 
the discrete speeds ci are nothing but the Gaussian knots, 
showing that Gaussian integration achieves a sort of au- 
tomatic "importance sampling" of velocity space which 
allows to capture the complexities of hydrodynamic flows 
by means of only a handful of discrete speeds. The LBE 
proves a very competitive tool for the numerical studies 
of hydrodynamic flows, ranging from complex flows in 
porous media to fully developed turbulence. 



B. Modified Lax-Wendroff scheme for species 
transport 

Since species transport equation is linear in the species 
concentration, we can solve it on a simple 6-neighbors 
cubic lattice. Within this approach, each species is as- 
sociated with a species density Cs, which splits into six 
separate contributions along the lattice links. 
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With these preparations, the transport operator in 3 
dimensions reads as follows (in units of = 1)): 

6 

Cs{x, t) = J^Pji^ - t - ^)Cs{x- Cj,t - 1) (12) 

3=0 

The index j runs over x and its nearest-neighbors (hence 
simpler than the LBE stencil) spanned by the vectors 
x + Cj,j = l,Q, j=0 being associated with the node x 
itself. The break-up coefficient pj represents the proba- 
bility that a particle at Xj = x — Cj at time t — 1 moves 
along link j to contribute to Cs {x) at time t. For instance 
in a one dimensional lattice the exact expression of these 
coefficients (in lattice units Cj = ±1, j = 1,2, At = 1) 
is: 

ft(x±i,i-i) = ^— + 2?;, i = i,2 (13) 

Po{x,t-l) = -2D'^ (14) 

where u' = (u + d^p) is the effective speed, in- 
clusive of the density gradient component, and D'^ = 
£>s(l — 'u'^)/2 is the effective diffusion, the square u' de- 
pendence being dictated by arguments of numerical sta- 
bility. 

C. Multiscale considerations 

The simulation of a reactive flow system is to all effects 
a multi-physics problem involving four distinct physical 
processes: 

1. Fluid Motion (F) 

2. Species Transport (T) 

3. Fluid- Wall interaction (W) 
4-. Wall Chemical Reactions (C) 

Each of these processes is characterized by its own 
timescale which may differ considerably from process to 
process depending on the local thermodynamic condi- 
tions. Loosely speaking, we think of F and T as to 
macroscopic phenomena, and W and C as of microscopic 
ones. The relevant fluid scales are the advectivc and 
momentum-diffusive time, and the mass-diffusion time 
of the species respectively: 




where L, H are the length and height of the fluid domain. 
The relevant time scales for species dynamics are: 

TD = H^/D, 

= K-\ (16) 
To = K-^ 



As discussed in the introduction, they define the major 
dimensionless parameters 

Re = UH/v = ta/t,, (17) 
Pe = UH/D = taItd. (18) 
Doc = Tc/ta, Doyj = Tw/ta (19) 

To acknowledge the multiscale nature in time of the 
problem, a subcycled time-stepper is adopted. This is 
organized as follows. The code ticks with the hopping 
time of the fluid populations from a lattice site to its 
neighbors dt = dx/c = 1. Under all circumstances dt 
is much smaller than both diffusive and advectivc fluid 
scales in order to provide a faithful description of fluid 
flow. Whenever dt exceeds the chemical time-scales (high 
Damkohlcr regime), fractional time- stepping, i.e. subcy- 
cling of the microscopic mechanisms, namely chemical- 
wall transfer is performed. This means that the chemi- 
cal and wall transfer operators are performed dt / Tc,dt / t^^ 
times respectively at each fluid cycle. As it will be ap- 
preciated shortly, since the flow solver ticks at the sound 
speed, the present microflow simulations proceed in very 
short time steps, of the order of tens of nanoseconds. This 
means that they can be in principle coupled to meso- 
scopic methods, such as kinetic Monte Carlo, affording 
a more realistic description of the fluid-wall interactions. 
In particular, a Kinetic Monte Carlo update of a sin- 
gle boundary cell could proceed in parallel with a cor- 
responding hydrodynamic treatment of the entire pile of 
fluid cells on top of the wall. The flip side of the medal 
is that in order to draw quantitative conclusions at the 
scale of the macroscopic devices a two-three decade ex- 
trapolation is required. This commands a robust scaling 
argument. 

IV. CATALYTIC EFFICIENCY: QUALITATIVE 
ANALYSIS 

Ideally, we would like to synthetize a universal func- 
tional dependence of the catalytic efficiency as a function 
of the relevant dimensionless numbers and geometrical 
design parameters: 

r] = f{Re,Da,Pe;g). (20) 

where g represents a vector of geometric parameters char- 
acterizing the boundary shape. The question is to as- 
sess the sensitivity of r; to ^ and possibly find an op- 
timal solution (maximum r]) within the given parame- 
ter space. Mathematically, this is a complex non-linear 
functional optimization problem for the geometrical pa- 
rameters. We flnd it convenient to start from a simple- 
and yet rcprcscntativc-bascline geometry as an "unper- 
turbed" zero order approximation, which is easily acces- 
sible either analytically or numerically. Perturbations to 
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this baseline situation can then be parametrized as "topo- 
logical excitations" on top of the geometrical "ground 
state". In the present study, the unperturbed geometry 
is a straight channel of size L along the flow direction and 
H xH across it. Perturbations are then defined as micro- 
corrugations in the bottom wall of the form z = h{x,y), 
h = being the smooth-wall unperturbed case. In this 
work, the perturbation is taken in the form of delta-like 
protrusions (barriers) h{x,y, z) = J2i hi5{x — Xi). 

From a macroscopic point of view the device efficiency 
is defined as amount of pollutant consumpted per unit 
mass injected: 

*m - ^out 



V 



where 



$(a;) 



j [uC]{x,y, 



z)dydz 



(21) 



(22) 



is the longitudinal mass flow of the pollutant at section x. 

The in-out longitudinal flow deficit is of course equal to 
the amount of pollutant absorbed at the catalytic wall, 
namely the normal-to-wall mass flow rate: 



r = / ^{x,y,z) ■ dS 
Is 



(23) 



where the flux consists of both advective and diffusive 
components: 



7 = uC - DVC 



(24) 



and the integral runs over the entire absorbing surface S 
The goal of the optimization problem is to maximize 
r at a given <I>i„. As it is apparent from the above ex- 
pressions, this means maximizing complex configuration- 
dependent quantities, such as the wall distribution of the 
pollutant and its normal-to-wall gradient. For future pur- 
poses, we find it convenient to recast the catalytic effi- 
ciency as 7/ = 1— T, where T is the channell transmittance 



(25) 



Prom a microscopic viewpoint, T can be regarded as 
the probability for a tracer molecule injected at the inlet 
to exit the channel without being absorbed by the wall 
and consequently it fixes the escape rate from the chemi- 
cal trap. Roughly speaking, in the limit of fast-chemistry, 
this is controlled by the ratio of advection to difi'usion 
timescales. More precisely, the escape rate is high if the 
cross-channel distance! walked by a tracer molecule in a 
transit time ta is much smaller than the channel cross- 
length H/2. Mathematically: Dta < i?^/4, which is: 



Pe > 4 L/H 



(26) 



The above inequality (in reverse) shows that in order 
to achieve high conversion efficiencies, the longitudinal 
aspect ratio L/H of the device has to scale linearly with 
the Peclet number. 



A. The role of micro-irregularities 

We now discuss the main qualitative eff'ect of geomet- 
rical roughness on the above picture from a microscopic 
point of view, i.e. trying to resolve flow features at the 
same scale of the micro-irregularity. 

In the flrst place, geometric irregularities provide a po- 
tential enhancement of reactivity via the sheer increase 
of the surface/volume ratio. Of course, how much of this 
potential is actually realized depends on the resulting 
flow configuration. 

Here, the fluid plays a two-faced role. First, geomet- 
rical restrictions lead to local fluid acceleration, hence 
less time for the pollutant molecules to migrate from 
the bulk to the wall before being convected away by 
the mainstream flow. This effect, usually negligible for 
macroscopic flows, may become appreciable for micro- 
flows with h/H ~ 0.1 (like in actual catalytic convert- 
ers), h being the typical geometrical micro-scale of the 
wall corrugations. Moreover, obstacles shield away part 
of the active surface (wake of the obstacle) where the fluid 
circulates at much reduced rates (stagnation) so that less 
pollutant is fed into the active surface. The size of the 
shielded region is proportional to the Reynolds number 
of the flow. On the other hand, if by some mechanism the 
flow proves capable of feeding the shielded region, then 
efficient absorption is restored simply because the pollu- 
tant is confined by recirculating patterns and has almost 
infinite time to react without being convected away. The 
ordinary mechanism to feed the wall is molecular diffu- 
sion/dispersion, which is usually rather slow as compared 
to advection. More efficient is the case where the fiow de- 
velops local micro-turbulence which may increase bulk- 
to- wall diffusive transport via enhanced density gradients 
and attendant density jumps — C^: 



^tur 



(27) 



where w' is the normal-to-wall microturbulent velocity 
fluctuation. This latter can even dominate the picture 

whenever turbulent fluctuations are sufficiently energetic, 
a condition met when the micro-Peclet number exceeds 
unity: 

Pe, = ^ » 1 (28) 

where h is the typical geometrical micro-scale. Given 
this complex competition of efficiency-promoting and 
efficiency-degrading interweaved effects it is clear that 
assessing which type of micro-irregularities can promote 
better eflSciency is a non-trivial task. 



B. EfBciency: analytic and scaling considerations 

For a smooth channel, the steady state solution of 
the longitudinal concentration field away from the inlet 
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boundary factors into the product of three independent 
one-dimensional functions: C{x,y,z) = X {x)Y (y) Z [z) . 
Replacing this ansatz into the steady-state version of the 
equation we obtain: 



X{x) = Xqc-^/' 

Y{v) = Fo 

Z{z) — ZoCos(z//_l) 



(29) 



with the longitudinal and cross-flow absorption lengths 
related via: 



/ = I' 



U 
D 



(30) 



where U is the average flow speed 



U{x) u{x, y, z)C{x, y, z)/ ^ C{x, y, z) (31) 



Note that the profile along the spanwise coordinate y 
remains almost flat because we stipulate that only the 
top and bottom walls host catalytic reactions. 

To determine the cross-flow absorption length l± we 
impose that along all fluid cells in a contact with the 
wall, the diffusive flux is exactly equal to fluid-to-wall 
outflow, namely: 



C _Cg 2 



T 



(32) 



where t the effective absorption/reaction time scale 

1 



1 1 

- ~ + 7 

T TD (Tc + Tu 



(33) 



and Nz = is the number of cells [dx — 1 in the code) 
in a cross-section x = const, of the channel. Therefore 
the factor 2/Nz is the fraction of reactive cells along any 
given cross-section x = const, of the channel. 

The form factor Cg/C is readily obtained by the third 
of Eq. (PI) which yields 



^ ~ cos(i7/2/j 



(34) 



Combining this equation with Eq. ( p^ ) we obtain a non- 
linear algebraic equation for Zj^: 



cos(A/2) 



2 



(35) 



where we have set A = H/l±. For each set of parameters 
this equation can be easily solved numerically to deliver 
l±, hence I via the Eq. (|3C|). 

Given the exponential dependence along the stream- 
wise coordinate x, the efficiency can then be estimated 
as: 



1 



(36) 



Note that in the low absorption limit L <^ I, the above 
relation reduces to 770 ~ L/l, meaning that doubling, 
say, the absorption length implies same efficiency with a 
twice shorter catalyzer. In the opposite high-absorption 
limit, L ^ the relative pay-off becomes increasingly 
less significant. 



C. Corrugated channel: Analytical estimates 

Having discussed tha baseline geometry, we now turn 
to the case of a "perturbed" geometry. Let us begin by 
considering a single barrier of height h. The reference 
situation is a smooth channel at high Damkohler with 
rjo = 1 — e^^/'. We seek perturbative corrections in the 
smallness parameter g = h/H, the coupling-strength to 
geometrical perturbations. The unperturbed wall- flux is 



2D^LH 
h 



(37) 



where Ch is the concentration at the tip of the barrier 
calculated in the smooth channel. Therefore Ch/h is an 
estimate of the normal-to-wall diffusive gradient. The 
geometrical gain due to extra-active wall surface is 



where 



El ~ ChUhhH 



Uh ^ 4C/o(5 - 5^) 



(38) 



(39) 



is the average longitudinal flow speed in front of the bar- 
rier along a section x = const.. The shadowed region of 
size w in the wake of the obstacle yields a contribution 



~ a D-^wH 



(40) 



where a is a measure of the absorption activity in the 
shielded region. 

Three distinctive cases can be identified: 

• a = 0: The wake region is totally deactivated, ab- 
sorption zero. 

• a = 1: The wake absorption is exactly the same as 
for unperturbed flow 

• a > 1: The wake absorption is higher than with 
unperturbed flow (back-flowing micro-vortices can 
hit the rear side of the barrier) 

Combining these expressions we obtain the following 
compact expression: 

Sr^ Fi+F2-F2(fe^0) Ah 

— = r - -77-FTR^h[Sc + K (a - 1)] 

770 i 2 H 

(41) 
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where A = H/L is the aspect ratio of the channel and 
Sc = v/D is the Schmidt number (fluid viscosity/tracer 
mass diffusivity) and the wake length can be estimated 
as w/h = KRch with K ~ 0.1. 

The above expression shows a perturbative (quadratic) 
correction in h over the unperturbed (smooth chan- 
nel situation). However, since the effective absorption 
in the shielded region is affected by higher order com- 
plex phenomena, the factor a may itself exhibit a non- 
perturbative dependence on h, so that departures from 
this quadratic scaling should not come as a surprise. 
Apart from its actual accuracy, we believe expressions 
like (|4^) may provide a qualitative guideline to esti- 
mate the efficiency of generic/random obstacle distribu- 
tions [xi,hi]: In particular, they should offer a semi- 
quantitative insights into non-perturbative effects due 
to non-linear fluid interactions triggered by geometrical 
micro- irregularities . 



the fluid- wall mass exchange is modelled via a mass trans- 
fer rate equation of the form previously discussed. 

We explore the effects of a sub-millimeter pair of barri- 
ers of height h a distance s apart on the bottom wall. The 
idea is to assess the effects of the interbarrier height, h, 
and interbarrier separation s on the chemical efficiency. 
Upon using a 80 x 40 x 40 computational grid, we obtain 
a lattice with dx — dy = dz = 0.0025 (25 microns), and 
dt — Cs dx/Vs ~ 50 10"^ (50 nanoseconds). Here we have 
assumed a sound speed Vs = 300 m/s and used the fact 
that the sound speed is = 1/ -s/S in lattice units. Our 
simulations refer to the following values (in lattice units): 
Uo ~ 0.1 - 0.2, D = 0.1, ly = 0.01, = = 0.1. This 
corresponds to a diffusion-limited scenario: 

= = 10 < TA ~ 800 <TD^ 16000 < = 160000 

(47) 

or, in terms of dimensionless numbers: 



V. APPLICATION: REACTIVE FLOW OVER A 
MICROBARRIER 

The previous computational scheme has been applied 
to a fluid flowing in a millimeter-sized box of of size 2 x 
1x1 millimeters along the x,y, z directions with a pair 
of perpendicular barriers of height h a distance s apart 
on the bottom wall (see Fig. |^ for a rapid sketch). 

The single-barrier set up corresponds to the limit 
s = 0. The fluid flow carries a passive pollutant, say 
an exhaust gas flow, which is absorbed at the channel 
walls where it disappears due to heterogeneus catalysis. 
The flow is forced with a constant volumetric force which 
mimics the effects of a pressure gradient. The exhaust 
gas is continuously injected at the inlet, x = 0, with a 
flat profile across the channel and, upon diffusing across 
the fiow, it reaches solid walls where it gets trapped and 
subsequently reacts according to a first order catalytic 
reaction: 



Pe ~ 40, Re ~ 400, Da > 80 



(48) 



C + A^P 



(42) 



where A denote an active catalyzer and P the reaction 
products. 

The initial conditions are: 



C(a;, y,z) ~ 1, x — 1 

C{x,y,z) ~ 0, elsewhere 
p{x,y,z) = 1 

u{x,y,z)^Uo, v{x,y,z) ^ w{x,y,z) 



= 



(43) 
(44) 
(45) 
(46) 



The pollutant is continuously injected at the inlet and 
released at the open outlet, while flow periodicity is im- 
posed at the inlet/outlet boundaries. On the upper and 
lower walls, the flow speed is forced to vanish, whereas 



As per the interbarrier separation, we consider the fol- 
lowing values: h/H = 0.2 and s/L = 0, 1/8, 1/4, 1/2, and 
h/H = 0.05,0.1,0.2 for s/L = 0. For the sake of com- 
parison, the case of a smooth wall (s — Q,h — 0) is also 
included. 

The typical simulation time-span is t = 32000 time- 
steps, namely about 1.6 milliseconds in physical time, 
corresponding to two mass diffusion times across the 
channel. The physico-chemical parameters given above 
are not intended to match any specific experimental con- 
dition, but rather to develop a generic intuition for the 
interplay of the various processes in action under the fast 
chemistry assumption. 



A. Single barrier: effects of barrier heigth 

We consider a single barrier of height h placed in the 
middle of the bottom wall at a; = L/2,z = 0. With 
the above parameters we may estimate the reference effi- 
ciency for the case of smooth channel flow. With C/ ~ 0. 1, 
and T — 20, we obtain I ~ 200, hence rjo ~ 0.5. 

A typical two-dimensional cut of the flow pattern and 
pollutant spatial distribution in the section y = H/2 is 
shown in Figs. || and ||, which refer to the case h = 8,s = 
(h/H = OA, s/L = 0.0). An extended (if feeble) recir- 
culation pattern is well visible past the barrier. Also, en- 
hanced concentration gradients in correspondence of the 
tip of the barrier is easily recognized from Fig. |^. A more 
quantitative information is conveyed by Fig. ^ where the 
integrated longitudinal concentration of the pollutant: 



C(x) = ^C(a:,y,z) 



(49) 
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is presented for the cases h = 0, 2, 4, 8 (always with 
s = 0). The main highhght is a substantial reduction 
of the pollutant concentration with increasing barrier 
height. This is qualitatively very plausible since the bulk 
flow is richer in pollutant and consequently the tip of the 
barrier "eats up" more pollutant than the lower region. 
In order to gain a semi-quantitative estimate of the chem- 
ical efficiency, we measure the the pollutant longitudinal 
mass flow: 



(50) 



The values at x ~ I and x ^ L define the efficiency ac- 
cording to Eq. ( ^ ) (to minimize finite-size effects actual 
measurements are taken at a; = 2 and x = 70). 

The corresponding results are shown in Table Jl where 
subscript A refers to the analytical expression (^) with 
a — 1. These results are in a reasonable agreement with 
the analytical estimate Eq. ( ^l| ) taken at a = f (same 
absorption as the smooth channel). However, for h — 8 
the assumption a — 1 overestimates the actual efficiency, 
indicating that the shielded region absorbs significantly 
less pollutant than in the smooth-channel scenario. In- 
deed, inspection of the transversal concentration profiles 
(Fig. H) along the chord x — 3L/4, y = H/2 reveals a 
neat depletion of the pollutant in the wake region. This 
is the shielding effect of the barrier. 

Besides this efficiency-degrading effect, the barrier 
also promotes a potentially beneficial flow recirculation, 
which is well visible in Figs. |^ and 0. Figure |^ shows 
the time evolution of the streamwise velocity u{z) in the 
mid-line x = 3L/4, y = H/2. It clearly reveals that recir- 
culating backflow only sets in for ft, = 8, and also shows 
that the velocity profile gets very close to steady state. 
A blow-up of the recirculating pattern in the near-wall 
back-barrier region is shown in Fig. ^ However these 
recirculation effects are feeble (the intensity of the recir- 
culating flow is less than ten percent of the bulk flow) 
and depletion remains the dominant mechanism. In fact 
for ft = 8 the measured local Peclet number is of the or- 
der 0.01 • 8/0.1 — 0.8, seemingly too small to promote 
appreciable micro-turbulent effects. In passing, it should 
be noticed that raising the barrier height has an appre- 
ciable impact on the bulk flow as well, which displays 
some twenty percent reduction due to mechanical losses 
on the barrier. 

Finally, we observe that the measured efficiency is 
smaller than the theoretical rjc for smoth channel. This 
is due to the fact that the flow ^{x — 2) is significantly 
enhanced by the imposed inlet flat profile C{z) = 1 at 
a; = 1 (as well visible in Fig. Leaving aside the initial 
portion of the channel, our numerical data are pretty well 
fitted by an exponential with absorption length I — 160, 
in a reasonable agreement with the theoretical estimate 
I ~ 200 obtained by solving Eqs. and dH). 



B. Effects of barrier separation 

Next we examine the effect of interbarrier separation. 
To this purpose, three separations s = 10, 20, 40 sym- 
metric around xq = L/2 are been considered. A typical 
two-barrier flow pattern with s = 40 is shown in Fig. ||. 
From this picture we see that even with the largest sepa- 
ration s = 40, the second barrier is still marginally in the 
wake of the first one. As a result, we expect it to suffer 
seriously from the aforementioned depletion effected pro- 
duced by the first barrier. This expectation is indeed con- 
firmed by the results reported in Table |l[ These results 
show that, at least on the microscopic scale, the presence 
of a second barrier does not seem to make any significant 
difference, regardless of its separation from the first one. 
As anticipated, the most intuitive explanation is again 
shadowing: the first barrier gets much more "food" than 
the second one, which is left with much less pollutant due 
to the depletion effect induced by the first one. Inspec- 
tion of the longitudinal pollutant concentration (Fig. |^) 
clearly shows that the first barrier, regardless of its lo- 
cation, "eats up" most of the pollutant (deficit with re- 
spect to the upper-lying smooth-channel curve is almost 
unchanged on top of the second barrier). Of course, this 
destructive interference is expected to go away for "well- 
separated" barriers with s ^ w. Indeed, the ultimate 
goal of such investigations should be to devise geomet- 
rical set-ups leading to constructive interference. This 
would require much larger and longer simulations which 
are beyond the scope of the present work. 



C. Effects of barrier height on a longer timescale 

Since the previous simulations only cover a fraction of 
the global momentum diffusion time, one may wonder 
how would the picture change by going to longer time 
scales of the order of /v. Longer single-barrier simula- 
tions, with t = 160, 000, up to 10 diffusion times, namely 
about 15 milliseconds, provide the results exposed in Ta- 



ble HI 



We observe that the quantitative change is very minor, 
just a small efficiency reduction due to a slightly higher 
flow speed. Indeed, the spatial distribution of the pollu- 
tant does not show any significant changes as compared 
to the shorter simulations, and a similar conclusion ap- 
plies to the flow pattern (see Figs. ^ and ^l|). This is 
because in a Poiseuille flow, the fluid gets quickly to, say, 
90 percent of its total bulk speed (and even quicker to 
its near- wall steady configuration), while it takes much 
longer to attain the remaining ten percent. Since it is 
the near-wall flow configuration which matters mostly in 
terms of a semi-quantitative estimate of the chemical effi- 
ciency, we may conclude that the simulation span can be 
contained to within a fraction of the global momentum 
equilibration time. 
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VI. UPSCALING TO MACROSCOPIC DEVICES 

It is important to realize that even tiny improvements 
on the microscopic scale can result in pretty sizeable cu- 
mulative effects on the macroscopic scale of the real de- 
vices, say 10 centimeters. Assuming for a while the ef- 
ficiency of an array of N serial micro-channels can be 
estimated simply as 

77^v - 1 - , (51) 

it is readily recognized that even low single-channel ef- 
ficiencies can result in significant efficiencies of macro- 
scopic devices with = 10 — 100 (see Fig. |l^). In 
particular, single-channel transmittances as high as 90 
percent can lead to appreciable macroscopic efficiencies, 
around 60 percent, when just ten such micro-channels are 
linked-up together. Such a sensitive dependence implies 
that extrapolation to the macroscopic scales, even when 
successfull in matching experimental data |ll|Jl^ , must 
be taken cautiously. In fact, the above expression (|5l| ) 
represents of course a rather bold upscaling assumption. 
As a partial supporting argument, we note that unless the 
geometry itself is made self-affine (fractal walls ||l^ ) , or 
the flow develops its own intrinsic scaling structure (fully 
developed turbulence), the basic phenomena should re- 
main controlled by a single scale I, independent of the 
device size L. Since both instances can be excluded for 
the present work, extrapolation to macroscopic scales is 
indeed conceivable. Nonetheless, it is clear a tight sin- 
ergy between computer simulation and adequate analyt- 
ical scaling theories is in great demand to make sensible 
predictions at the macroscopic scale. 
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VII. CONCLUSIONS 



This work presents a very preliminary exploratory 
study of the complex hydro-chemical phenomena which 
control the effective reactivity of catalytic devices of mil- 
limetric size. Although the simulations generally confirm 
qualitative expectations on the overall dependence on the 
major physical parameters, they also highlight the exis- 
tence of non-perturbative effects, such as the onset of 
micro-vorticity in the wake of geometrical obstrusions, 
which are hardly amenable to analytical treatment. It is 
hoped that the flexibility of the present computer tool, as 
combined with semi-analytical theories, can be of signifl- 
cant help in developing semi-quantitative intuition about 
the subtle and fascinating interplay between geometry, 
chemistry, diffusion and hydrodynamics in the design of 
chemical traps, catalytic converters and other related de- 
vices. 
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Run 


h/H 


n 


t>V "VA 
V ' VA 


ROO 





0.295 


0.00 


R02 


1/20 


0.301 


0.02,0.025 


R04 


1/10 


0.312 


0.06,0.10 


R08 


2/10 


0.360 


0.22,0.40 



TABLE I. Single barrier at a; = 40: the effect of barrier 
height. 



Run 


s/L 


V 


ROO 





0.30 


R08 


1/8 


0.36 


R28 


2/8 


0.37 


R48 


4/8 


0.375 



TABLE II. Two barriers of height h = 8: Effect of inter- 
separation s. 



Run 


h/H 


'/ 


ri • r,4 


LOO 





I). 290 


(),() 


L02 


0/20 


0.296 


0.02,0.025 


L04 


1/10 


0.307 


0.06,0.10 


LOS 


2/10 


0.360 


0.24,0.40 



TABLE III. s = 0, /i = 0, 4, 8: 10 mass diffusion times 
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FIG. 1. Sketch of the of a section at y = const, of a typical 
channel with two microbarriers. Two barriers of height h = 3 
a distance s = 10 apart: F=fluid, B=buffer. 

u(x,z) at y=U2: t=32000 
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FIG. 2. Typical two-dimensional cut of the flow pattern 
with a single barrier of heigth /i = 8. Streamwise flow speed 
in the plane y = H/2. 
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C(x,z) at 7=172: t=32000 
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FIG. 3. Concentration isocontours with a single barrier of 
heigth h = 8. 



Transverse pollutant concentration at x=3L/4,y=H/2: 32000 steps 
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FIG. 5. Transverse pollutant concentration C(z) at 
X = 3L/4 and y = H/2. Single barrier of varying height. 
The four curves for each of the three different heigths are 
taken at t = 3200, 6400, 29800, 32000. 



Longitudinal pollutant concentration: single tiarrier 
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FIG. 4. Integrated longitudinal concentration C{x) of the 
pollutant with a single barrier of height h = 8 after 32000 
steps. 



streamwise speed U(Z) at x=3i-/4,y=L/2 
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FIG. 6. Time evolution of the transversal streamwise speed 
u{z) at a; = 3L/4 and y = L/2. Single barrier of varying 
height. 
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stream function at y-L/2: h=8, 160000 timesteps 



FIG. 9. Longitudinal concentration C{x) for h = 8, 
s = 0, 20, 40 all at t = 32000. 
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FIG. 7. Blow-up of the streamlines of the flow field past 
a barrier of height h = 8 located at x = 40. The velocity 
direction in the closed streamlines of the vortex is clockwise. 



Streamwise u(x,z) at y=L/2 



Transversal pollutant concentration at x=3L72,y=H/4: 1 60000 time-steps 
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FIG. 10. Integrated longitudinal concentration C{x) of the 
pollutant with a single barrier of height h = 8 after 160000 
steps. 
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FIG. 12. Efficiency of a scries of micro-channels as a 
function of the single-channel transmlttance. 



